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We present a complete-quantum description of multi-particle Szilard engine which consists of a 
working substance and a Maxwell's demon. The demon is modeled as a multi-level quantum system 
with specific quantum control and the working substance consists of identical particles obeying 
Bose-Einstein or Fermi-Dirac statistics. In this description, a reversible scheme to erase the demon's 
memory by a lower temperature heat bath is used. We demonstrate that (1) the quantum control 
of the demon can be optimized for single-particle Szilard engine so that the efficiency of the demon- 
assisted thermodynamic cycle could reach the Carnot cycle's efficiency; (2) the low-temperature 
behavior of the working substance is very sensitive to the quantum statistics of the particles and 
the insertion position of the partition. 



I. INTRODUCTION 

Maxwel introduced, in 1871, a notorious being, known 
as Maxwell's demon nowadays, to discuss the "limita- 
tions of the second law of thermodynamics" [l[ . Such 
a demon distinguishes the velocities of the gas particles 
and controls a tiny door on a partition of the gas con- 
tainer to create a temperature difference, which breaks 
the Clausius statement of the Second Law of Thermo- 
dynamics (SLoT). To reveal the essence of the Maxwell's 
demon, Leo Szilard proposed a single-particle heat engine 
0, named as Szilard heat engine (SHE). The demon in 
SHE distinguishes the positions of the particles after the 
partition has been inserted. With the help of the de- 
mon, SHE can absorb heat from a single heat resource 
and convert it into work without apperant evoking other 
changes in the cycle. Szilard pointed out that the SLoT 
was no longer violated if one considered the entropy in- 
crease during the measurement. Brillouin generalized the 
Szilard's argument and identified the thermodynamic en- 
tropy with the informational entropy firstly [3| . 

However, the measurement could be carried out with- 
out any change in entropy 4|. Acutally, it was realized 
that a logically irreversible process must be "accompa- 
nied by dissipative effects" [5[ in the physical realization 
of the information processing, which is known as Lan- 
dauer's erasure principle. Bennett used this point of view 
in the study of the Maxwell's demon paradox and pointed 
out that the erasure of the demon's memory instead of 
the measurement was logically irreversible and thus must 
be accompanied by dissipative effects [il] . With these ob- 
servations, the conventional cycle presented by Szilard is 
indeed not a thermodynamic cycle because the demon's 
memory has not been erased to complete the cycle. The 
SLoT will be saved if one considers the erasure process 
to finish the cycle of the demon. As people believe the 
essence of information should be discussed in the frame- 



* Elec tronic address: |suncp@itp.ac.cn| 

URL: |http: //power ■ itp. ac . cn/suiicp/lndex.html| 



work of quantum mechanics, various quantum versions 
of SHE have been proposed with different views about 
quantum measurements. One proposal is semi-classical 
[l3| • The working substance in this proposal is quan- 
tum mechanic while the demon is considered as a classical 
controller whose role is to extract information through 
measurement and control the system. The paradox of 
Maxwell's demon was solved by arguing Landauer's era- 
sure principle. However, the final solution should include 
the MD in the cycle and treat also the MD in a quantum 
fasion0, [in. For SHE, it is also proved the existence of 
MD will not violate the SLoT in Ref. [8], where MD is 
modeled as a two-level system. 

The focus of the study in both classical and quantum 
mechanical frameworks is the erasure process, which is 
crucial to solving the Maxwell's demon paradox. In an 
odinary way, the demon and the working substance are 
in contact with the same heat bath. After erasing the 
demon by applying some work, one will find that SHE 
can not extract work in a cycle at all and the SLoT is 
not violated. However, a more general erasure should be 
done with a lower-temperature heat bath which is also 
called a heat sink. In this suituation, we turn the SHE 
into a thermal dynamic cycle, where the non-violation 
of SLoT can be proved by illustrating non-exceeding of 
Carnot's efficiency. It is realized that the effective tem- 
perature of the MD's initial state actually charasticrizes 
the error in the control of the heat engine0, Our 
previous work in Ref. iS] emphasized the functions of 
the demon with errors in the study of the siiigle-particle 
SHE. But the erasure schemes in Refs. 0, Q are irre- 
versible. Therefore, the efficiencies of the heat engines in 
these papers can not reach the Carnot cycle's efficiency. 
One purpose of the present paper is to establish an opti- 
mal scheme of the thermodynamic cycle with a reversible 
erasure process, assisted by demon. It is shown that the 
partition-removing process is not always reversible, which 
leads to the lower efficiency. We aso find the existence 
of the optimal expansion position to improve the effi- 
ciency of the single particle SHE to the Carnot cycle's 
efficiency. The other purpose of this paper is to reveal 
the role of the quantum statistical properties of the work- 
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ing substance. We generalize our previous works about 
demon-assisted quantum heat engine by using a multi- 
particle working substance, which is the ideal Bose or 
Fermi gases, and find that the ratio of the work extracted 
to the working temperature has some discontinuous be- 
havior, and that discontinuous behavior is closely related 
to the degenerate-ground-state phenomenon. 

The paper is organized as follows: In Sec. |IT] we de- 
scribe the model of quantum multi-particle SHE and the 
working scheme briefly. In Sec. IIIIl we study in details 
the five steps of the working scheme: insertion, measure- 
ment, controlled expansion, removing and erasure sep- 
arately and calculate the work applied and heat trans- 
ferred in each step. In Sec. lIVl the efficiency of the engine 
is evaluated. It is found that our heat engine's efficiency 
can not exceed the Carnot cycle's efficiency. For single- 
particle SHE, we optimize the scheme to make the effi- 
ciency of the engine reach the one of Carnot's. In Sec.lVl 
we discuss the behavior of the engine in low-temperature 
regime and show the ratio of the work extrated to the 
working temperature is closely related to the degener- 
ate point (whose definition can be seen in this section) 
and the particles' statistical properties. Conclusions and 
remarks are given in Sec. IVIl 
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Figure 2: (Color Online) The states of the chamber and the 
demon (a) before and (b) after the expansion process, (a) 
After the measurement, the demon has revealed the number 
of the particles in the right compartment i and recorded this 
information in its memory. If the demon is error-free, then 
its inner state (or memory) will be However, when the 
demon contains some errors, there would be some mistakes in 
the demon's memory. For example, the sub-figure(2) shows 
a situation when z = and the demon is in state jl). Thus, 
there are total (A^ + 1)^ situations after the measurement for a 
system with particles. Here, the total 9 situations are listed 
for N — 2. (d) After the measurement, the partition acts as 
a piston under the control of the demon. If the demon is in 
the state the piston will move to a new position h. Such 
a movement is called controlled expansion. Here the results 
of the expansion for all 9 situations are listed for N = 2. 



Figure 1: (Color Online) Modeling the chamber of gas by 
a one-dimensional infinite square well, (a) The chamber is 
modeled as a one-dimensional infinite square well with width 
L where the particles are confined. The inner lines in the 
well represent the single-particle energy levels of the well, (b) 
After the insertion, the chamber is split into two chambers, 
which are modeled as two 1-D infinite square wells with widths 
I and L — I respectively. The inner lines represent the single- 
particle energy levels. 



II. QUANTUM MULTI-PARTICLE SZILARD 
ENGINE 

The working substance is modeled as a collection of N 
particles confined in a chamber, which is described as a 
one-dimensional infinitely deep square well with width L 
illustrated in Fig. [Ija). Generally, the particles can sat- 



isfy any distribution, Bose distribution or Fermi distribu- 
tion. We calculate the canonical partition function Zn{l) 
as follows, where the subscript n denotes the number of 
the particles in the well and I represents the width of 
the well. Its explicit expression is given in the Appendix 
El In the following calculations, we do not specify the 
concrete partition functions in most cases. The results 
obtained are valid for both Bose and Fermi distributions 
except for the case we claim specifically. 

In the model. Maxwell's demon is described as a quan- 
tum system with + 1 energy levels to record the com- 
plete information about how many particles in each com- 
partment after the insertion process. We may use a de- 
mon with number of energy levels less than A'^ -I- 1, but it 
can only record partial information of the working sub- 
stance, which decreases the work extract. We consider 
in the present paper only the demon with A^ -|- 1 levels, 
with \i) be its i-th eigen-state and be the energy of 
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state \i). The energy of |0) is set to be the zero point 
of energy, i.e., Aq — 0. For a demon without errors, its 
initial state is a pure state, namely, pdo = |0) (0|. In this 
paper, we study the demon with errors generally, whose 
initial state can be written as 

Pdo = P^o^ |0) (0| |1) (1| + • • ■+p^^^ \N) {N\ . (1) 

When pf''' vanish for all i 0, this demon returns to the 
one without errors, wherejO) is the standard state. For a 
demon with errors, a larger population pg"' in the state 
0) means a more efficient demon Q. 

The thermodynamic cycle consists the five steps. At 
the beginning, the working substance is in equilibrium 
with the heat resource at a higher temperature Ti, and 
the demon is initially in the state pdo (see Eq. ([T])). In 
the first step, a partition is inserted at the position /, 
which splits the initial well into two with widths being 
I and L — I respectively (see Fig. (Hb)). After the inser- 
tion, there are + 1 possible situations and one does 
not know which situation happens. Suppose i to be the 
number of the particles in the right compartment. Then 
the A^ + 1 situations correspond to i = 0, 1, ■ • • , A^. In 
the second step, the demon measures the working sub- 
stance and records this information in its memory. For 
the error-free demon, it would be in the state \i) for the i- 
th situation after the measurement. However, the demon 
with errors may record wrong information. For example, 
the working substance is in 0-th situation while the de- 
mon is in the state |1) as illustrated in Fig. ^a). In the 
third step, the partition acts as a piston and expands un- 
der the control of the demon. Specifically, if the demon 
is in the state \i), it will move the piston to the position 
k (see Fig. [^Jb)). Actually, one can optimize the expan- 
sion position li to get the best efficiency. In the fourth 
step, the partition is removed and the working substance 
returns to its initial state. After these four steps, the 
working substance absorbs heat from the heat resource 
and converts it into work. It seems to violate the SLoT. 
However, the demon has not returned to its initial state 
to complete the cycle. Thus, an additional step, step five, 
is needed to erase the memory of the demon. The effect 
of this step saves the SLoT as pointed out by Landauer's 
erasure principle. 



III. QUANTUM THERMODYNAMIC CYCLE 
OF MULTI-PARTICLE SZILARD ENGINE 

In this section, we analyze in detail the thermodynamic 
cycle of the quantum multi-particle SHE. We calculate 
the heat absorbed and work applied in each step. 



heat resource and the process is quasi-static. Before the 
insertion, the state of working substance pso is 



Pn{L) 



1 



Zn{L) 



-/3iffjv(i) 



(2) 



where /3i = l/Ti is the inverse temperature and the 
Boltzmann constant is set to be unit, i.e., /c^ = 1, and 
Hn{L) represents the Hamiltonian of A^ particles con- 
fined in the infinite square well with width L. Z^{L) is 
the corresponding partition function and its explicit ex- 
pression is presented in Appendix [Xj After the insertion, 
the chamber is split into two parts with widths I and L — I 
respectively and the substance is in an equilibrium state 
as 



N 

E 

i=0 



P.{l)p'^,^S)®pf{L^l). 



(3) 



Here, p'^_^{l) describes an equilibrium state of the left 
chamber with similar meaning as those in pn{L). 



m) 



ZN-^{l)Z,[L-l) 

J:1,Zn^,{1)Z,{L~1) 



(4) 



is the probability to find i particles in the right chamber 
with insertion position being I. As the process is isother- 
mal, the work applied in this step is the difference of the 
free energies after and before the insertion, i.e.. 



and the heat absorbed is 

Qins — ^l(5's,i] 



sO, 



Sso)- 



(5) 



(6) 



The free energy and the entropy are given in terms of the 
partition functions as 



-Tiln 



'ZN-^il)Z,iL-l) 



F,o = -TilnZAr(L), 



(7) 
(8) 



and 



^'^'^'Wi^^'' (^5]z^-.(0^.(i-0^ ,(9) 



Sso = (l-/3i^)lnZjv(i), 



(10) 



respectively. As shown in Ref . , in quantum mechanics 
framework, the insertion work Wins is no longer zero. 
Moreover, it was proved in Ref. Q that limT->oo Wins — 
oo. 



step 1: Insertion 



Step 2: Measurement 



In the first step, the partition is inserted isothermally, 
which means the working substance is in contact with a 



In the second step, the total system is isolated from 
the heat bath. The demon finds out the number of the 
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particles in the right compartment. This measurement 
process aims at estabhshing the correlation between the 
working substance and the demon. For the demon with 
no error, the state of the demon after the measurement 
is \i) when there are i particles in the right compart- 
ment. For the demon with errors, its initial state is a 
mixed state. Thus, we should appoint the final states 
of the demon for all situations. Let the final state be 
\fi{j)) when the demon is initially in the state \j) and 
the number of the particles in the right compartment is 
i. Mathematically, fi'-j^ fi{j) is a map from the set 
{0, 1, • • • , N} to itself. If the demon is error- free, the fi- 
nal state will be \i) when there are i particles in the right 
compartment, which leads the first constraint of /i, i.e., 
/i(0) — i. With this notation, the operator representing 
this measurement is 



N 



N 



(11) 



N N 



where ) represents the /^-th eigen-state of the working 
substance when there are i particles in the right compart- 
ment. For a physical operation, U should be unitary (so 
that it is a quantum nondemolition pre-measurement), 
which leads to another constraint of fi{j). That is, for a 
fixed i, when j runs over from to N, fi{j) should also 
run over from to N. Thus, {fi{j)\j = 0, 1, • • • ,N} is 
a permutation of {0, 1, • • ■ , N} which satisfies /i(0) = i. 
When = 1 , there is only one function of such kind and 
the corresponding operator U represents a controUed- 
NOT operation. For > 1, such function exists. One 
realization is fi{j) = i + j (mod iV + 1). 



After the measurement, the state of the total system 
becomes pmoa = C^Ps.ins ® PdoU\ namely 



3=0 



i=0 



I), 



(12) 



where : j i— ^ fi~^{j) is the inverse of the map fi. It 
follows from Eq. (|12l) that there is a possibility that the 
demon is in the state \i) while the number of particles in 
the right compartment is not i. That is exactly what we 
mean 'error' within the demon. Due to the isolation from 
the heat bath, there is no heat transferred, Qmea = 0, 
and the work applied is just the energy change during 
the measurement: 



IS 



N 



w=y: 



J=0 



(13) 



which compensates for the energy difference between the 
initial state and the final state of the demon during the 
measurement. Therefore, if all the energy levels of the 
demon are degenerate, one will not apply any work to 
perform the measurement. Thus the work applied may 
not always be non-zero 



Step 3: Controlled Expansion 

In this step, the partition in the chamber acts as a pis- 
ton and will be moved to the position according to the 
memory of the demon. Specifically, if the state of the de- 
mon is \i), the finial position will be k, which is called ex- 
pansion position. This expansion process is slow enough 
and the working substance is in contact with the heat 
resource during the total expansion process. Thus the 
state of the total system after the controlled expansion 



N 



N 



(14) 

However, the 'impenetrability' of the piston does not 
make this process an isothermal process since the tran- 
sition between two states with different numbers of the 
particles in the right compartment is forbidden. Thus we 
should deal with this expansion process for each situation 
separately. For the situation when the demon is in state 
\j) and there are i particles in the right compartment, 
the work applied Wij is 



-Ti In ZN-^ {lj)Z^ (L Ti \nZN-^{l)Z^{L-l), 



and the heat absorbed Qij is 

Qri = T^{s{p%^Sj)®pfiL-ij))-s{p%^dmpfiL-i^^^^ 

The probability of this situation is pij = Pi{l)p^ji^li^^y 

Thus, the total work applied and the total heat absorbed 
are: 



N N 
N N 

Qcxp — ^ ^ ^ ^ Pij Qij : 

j=0 i=0 



(15) 
(16) 
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respectively. One can prove that even if the controUed 
expansion is not isothermal, the heat absorbed also satis- 
fies the relationship Q = TAS", where AS* is the entropy 
change of the total system. Since the states with differ- 
ent numbers of the particles in the right compartment 
are orthogonal to each other, the entropies of the total 
system before and after the controlled expansion are 



Step 4: Removing 



(17) 



N N 

j=a 1=0 
N N 

j=a i=o 

respectively, where H{{pij}) = -J2ijPij^^Pv- Thus, 
we have the following relationship 

TiAS — Tl(S'oxp — S'lnca) 
N N 

j=0 i=0 



N N 



^ ^ ^^PijQij Qcxp- 



(19) 

















• • 

• • 


• 

• 

• 




• 

• • 
• • 
• 


• • 




• • • 

• • 

• 




(a) 
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(c) 



Figure 3: (Color Online) The removing process, (a) Before 
the removing, the state of the working substance may not be 
an equilibrium state because of the 'impenetrability' of the 
partition (see the sub-section (b) When the partition is 
lifted a little, it is no longer impenetrable and the particles 
can fly between the two sides of the chamber freely to reach 
the equilibrium state, (c) The partition keeps being lifted and 
finally be removed isothermally. 



After the controlled expansion, the partition will be 
removed slowly while the working substance is in contact 
with the heat resource. Because the state in Eq. 
before removing is not a thermal equilibrium state (it 
is because the controlled expansion is not an isothermal 
process, see the last subsection), one should carefully deal 
with this step. Once lifted a little, the partition is no 
longer 'impenetrable' and the particles can fly between 
the two sides of the chamber freely to reach an equilib- 
rium state. Thus the removing process consists of two 
sub-steps. After the first sub-step finished, the parti- 
tion has been lifted to create a small slit with width dp 
through which the particles can fly from one side to the 
other. Whatever small the width dp of the slit is, the 
working substance can reach its equilibrium state. Thus, 
we can let the width of the slit tend to zero, i.e., dp 0, 
which makes this sub-step a thermalization process of 
the working substance while the partition remains still. 
Thereafter the state of the system after the first sub-step 
becomes 



N N 

Prev - 1^') 01 ® E Pdh)PN-^ilJ} ® pf " Ij), (20) 

j=Q i=a 

where p^^^ = 'Yl!i=oPiQ')P f-'^(j)- During this process, there is no work applied (since dp 0) and the heat absorbed 
is exactly the difference between the inner energies of the total system before and after the first sub-step, 

N N 

Qrov = E E [pf^P^i^j) - [U{p^N-^m + U{pf{L - I,))] . (21) 

After the thermalization, the working substance is in a thermal equilibrium state, which makes the rest sub-step, the 
second sub-step, an isothermal process. After the removing, the working substance returns to its initial state and the 
correlation between the working substance and the demon no longer exists. During the second sub-step, the work 
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applied and the heat absorbed are 

Qicv — ^1 



N / N \ 

InZAr(L) + ^pf ^ In ZN-^{h)Z^{L - /,) 

j=0 \j=0 / 



N 



j=0 



(22) 

(23) 



respectively, and the state of the total system after the 
removing process is 



concrete, we assume that p'^^ = 0.7 while p^'' — 0.2 
and pj"'' — O-l- Suppose that the populations become 



.(0) 



N 



(24) 



One can see that the state of the demon does not return 
to its initial state, which does not make the above four 
steps a thermodynamic cycle. Thus, an addition step, 
erasure, is needed. 



Step 5: Erasure 

In Eq. ([M)) . the state of the demon after the removing 
process is 



N 



Pdl 



(25) 



which is different from the initial state of the demon. As 
part of the engine, the demon should be re-initialized to 
complete the cycle. Because the entropy of pdi is differ- 
ent from the entropy of pdo, such a re-initialization pro- 
cess, or erasure process, can not be realized by a unitary 
evolution and must be a logical irreversibility process ac- 
companied by dissipative effects A natural way to 
complete the re-initialization process is to make use of a 
thermalization process. It can be shown that the heat 
bath used in the thermalization process is a heat sink 
instead of a heat resource. Because the entropy of pdi is 
larger than the entropy of pdo, the heat must be trans- 
ferred from the demon to the heat bath. It seems that 
after a series of complex tasks, the demon needs to be 
cooled down. According to thermodynamics, the best re- 
initialization process is a reversible one which does not 
cost additional dissipative work. Before going into the 
general discussion, we first consider a three-level demon 
as an example to find out how to re-initialize or erase the 
demon reversibly. 

The three-level demon considered is shown in Fig. 21 
where the three states are |0), |1) and |2) respectively. 
The energy of the state |0) is set to be zero and the en- 
ergy of the state \i) is for z = 1, 2. For a demon in its 

initial state, the populations of the three levels are pg^"* , 
pf'' and p^2 ^ respectively (see Fig. IHJd)). To be more 



0.3, p'^i^ = 0.5 and P2' = 0.2 respectively after 
the removing (see Fig. Ufa)). In order to re-initialize it 
by using a reversible scheme by means of a heat bath 
with T2 = 1, we adjust the demon's energy level spacings 
adiabatically (see Fig. |4](a)-(b)). After being adjusted, 
the level spacings are changed from (i = 1,2) to A,^ 
(see Fig. Hfb)). During this process, the population of 
each state remains unchanged. Thus, the new level spac- 
ings satisfy the following relationship to ensure that the 
effective temperature of the demon is T2, i.e.. 



p11_ 



exp 



T2 



(26) 



for I = 1,2, which implies A^ = -0.5108 and A^ = 
0.4055. Next, the demon is kept in contact with the 
heat bath at temperature T2 without any undesired irre- 
versible heat transfer, and we adjust the energies of its 
levels isothermally (see Fig. |l](b)-(c)). In this process, 
the population of each level changes as the level spacings 
change and the final populations are adjusted to be p^^\ 

p^i^ and p^2'' respectively, which are the initial popula- 
tions of the demon. Therefore the new level spacings A" 
{i = 1,2) satisfy the relationship 



(0) 

Pi 

(0) 

Pi) 



exp 



T2 



(27) 



which implies A'/ = 1.2528 and A^' = 1.9459 (see Fig. 
Iljc)). After that, the population of each state has re- 
turned to its initial value while the level spacings do not. 
Thus we need another adiabatical adjusting process as 
in the first step to re-initialize the level spacings (see 
Fig. 2] (c)-(d)). These three sub-steps complete the re- 
initialization process and are all reversible, which makes 
the total re-initialization process reversible. 

The above example inspires us to construct the re- 
versible re-initialization scheme in a more general case. 
In the following discussion, a detailed scheme of the re- 
initialization process is presented. In the first sub-step, 
we adjust the energy level of the state |j) for j 7^ adia- 
batically while the energy of |0) is kept fixed. The corre- 
sponding new state is denoted as During this sub- 
step, the demon is isolated from the heat sink. The target 
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(d) 



A 



step3 



step2 




(b) 



Figure 4: (Color Online) Erasing the memory of a three- 
energy-level demon reversibly. Suppose that the populations 
of the three levels after the removing process the partition are 



Pq^^ — 0.3, p'l''' = 0.5 and pj^'' = 0.2 as in (a) and the initial 



state's populations of the three levels are Pq"' = 0.7, pf = 0.2 
andp^°' =0.1 as in (d). In step 1 (a-b), the demon is isolated 
from the heat bath and the energy level spacings are adjusted 
from Ai to for i = 1,2. During this process, the popula- 
tions of the levels remain constant. The new level spacings 
are chosen so that the demon has an effective temperature 
T2 — 1 after the adjusting, which implies A'^ — —0.5108 and 
A2 — 0.4055. In (b), the minus sign in front of A'l is due 
to the negative value of A'^. In step 2 (b-c), the demon is 
in contact with a real heat bath at temperature T2 and the 
energies of the levels are adjusted. In this step, the popula- 
tions of the levels change while its energy spacings change. 
The goal of this step is to let the populations of the levels 
return to be pg'^ , pj"' and pj"' . Thus the new level spacings 
are A'/ = 1.2528 and A^' = 1.9459. In step 3 (c-d), the demon 
is isolated from the heat bath again and its level spacings are 
adjusted to return to its initial values. The isolation guaran- 
tees that the populations of the levels are not changed. Thus 
the state of the demon has returned to its initial state. 



is to let the demon's effective temperature be T2, which 
is the temperature of the heat sink. Thus, the new level 
spacing satisfies 

exp(-/32Aj-) = ^ , (28) 

and the demon after the first sub-step is in the state 



N / N \ 

j=0 \i=0 ' / 



(29) 



It should be noticed that Eq. ((28)) always has a solu- 
tion except for some extreme conditions, for example, 
J2iLo ^i(^)Pf^^{j) ~ ^ ^2 = 0. If these extreme con- 
ditions happen, we will use some tricks to deal with it. 
In the case of T2 = 0, we will introduce a small positive 
quantity e and let be T2 + e. Next we will use 
instead of T2 to calculate the parameters of the scheme 
such as A- and A". Such a scheme is not reversible. 
However, if e tends to zero, the scheme will tend to be- 
have as a reversible scheme. In the following, we assume 
that Eq. ([28]) does have solution. In this sub-step, there 
is no heat transfer and the work applied is 




m)p 



fiij) 



(A^-A,). 



(30) 



In the second sub-step, demon is in contact with the 
heat sink and the j-th level's energy is adjusted from A^- 
to A" isothermally, where A" satisfies 



-/32A^' 



Pi 



(31) 



and is the inverse temperature of the heat sink. Thus, 
the state of the demon after this sub-step is 



N 
J=0 



(32) 



where represents the state corresponding to after 
the adjusting. The work applied and the heat transferred 
during this process are 



1 ^ 

— [ln(l + e-l^^^'^) - ln(l + e-'^^'^")] (33) 



Qini — T2 



respectively. Here H{{pi}) — — ^-pi Inpi represents the 
Shannon entropy of the probability distribution {pi}. 

In the last sub-step, the demon is isolated from the 
heat sink again and the energies of all of its levels are 
adjusted to its initial values. Thus the final state of the 
demon is J2jPj \j) (j Ii which is indeed its initial state. 
During this sub-step, no heat transfers and the work ap- 
plied is 



ini,3 



N 

E 



Pji^j 



(35) 



These three sub-steps together make up the re- 
initialization scheme. Such a scheme involves adjusting 
the energies of the demon's levels, which is not easily re- 
alized. In a concrete experiment, if the quantum system 
that functions as the Maxwell's demon has unchangeable 
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energy levels, we will only use an irreversible scheme as 
in Refs. 0, i- As we know in thermodynamics, such 
an irreversible re-initialization scheme reduces the effi- 
ciency of the heat engine because of the dissipative ef- 
fects. Thus, when it comes to the Maxwell's demon 
paradox, one needs a reversible re-initialization scheme 
without dissipative effects to get a good physical picture. 



3. The work extracted is 



W 



-Wins - 



- VFi„i.2 



cxp ^^rc 



(38) 



IV. 



EFFICIENCY OF SZILARD ENGINE 



Combining the results obtained in the above section, 
one finally gets the total heat transferred and the work 
extracted in the thermodynamic cycle: 

1. The heat absorbed from the heat resource is 

Ql — Qins + Qoxp + Qrcv + Qrcv (36) 

N N 

= Ti H({/',(0}) + ^^P,(OP/-io)lnP,(/,) 

2. The heat released to the heat sink is 

Q2 = -Qini (37) 



= T2 



N 



Hi{j2m)Pf-^u)^)~H{{p,}) 



1=0 



which can be checked to satisfy the first law of thermo- 
dynamics, i.e., W = Ql — Q2- 



Firstly, it is emphasized that the result is general with- 
out referring to the statistical properties of the working 
substance, since Pi{l) has contained all the information 
of the distribution. Thus Eq. §^,Eq. ^ and Eq. §Eii 
are correct both for the Bose system and the Fermi sys- 
tem. Secondly, one can optimize the expansion position 
Ij to maximize the efficiency of SHE. Consider a two- 
particle SHE as an example, in which the demon is of 
three-level. The function fi{j) we choose is fi{j) = i + j 
(mod 3) . Thus the heat absorbed from the heat resource 



Ql = TiiJ({P,(0}) +PQ{l)po\nPo{lQ) + Pi{l)p2 \nPi{h) + P2{l)pi lnP2(?o) 
+Poil)pi InPoih) + Piil)palnPi{h) + P2{l)p2 \nP2{h) 
+Po{l)p2 InPoih) + Pi{l)pi InPiih) + P2{1)PQ lnP2(/2). 

Thereafter, the optimized expansion position U should satisfy 

dQi dQi dQi 



dlo dli 



dh 



= 0. 



(39) 



(40) 



Working out the above equation, one will get the optimized expansion positions Zomax, ^imax and Z2max- 

Thirdly, an upper limit of the heat absorbed from the heat resource Qi in Eq. (j36|) is estimated. Use the method 
presented in the Appendix |BJ one gets 



Ql ^ Ti 



N N 



i=0 j=0 



^({E^^^P/r^O)}^ -H{{p^}) 



(41) 



Therefore the efficiency of the SHE has the upper limit: 

, Q2 

^1-^- (42) 



Thus, it is clear that the efficiency of the SHE can not 
exceed the Carnot cycle's efficiency when the erasure pro- 
cess or re-initialization process is considered, which saves 
the SLoT. It can be seen more clearly when = 1, which 
represents a single-particle SHE. For a single-particle 
SHE, there is only one realization of /i(j), and the cor- 
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responding measurement process is indeed a controUed- 
NOT operation as in Ref. 0,[1], i.e., /o(j) = j for j — 0,1 
while /i(l) = and /i(0) = 1. The formulas Eq. §Q 
and Eq. ([57)) are reduced to 



Qi=Ti 



H{Po{l),Pi{l)) 
+Po{l)po InPo(^o) + Pi{l)Pi lnPi(?o) 
+Po {l)pi \nPo{li) + Pi {l)po \nPi{h) 



and 



H{Po{l)pi+Pi{l)po,Pa{l)po + Pi{l)Pi) 
-H{p„,pi) 



respectively. It is easy to maximize Qi to find that lon 
and hmax satisfy 



max I 



Poil)po 



Po{l)Po + Pi{l)pi 

Po{l)pi 
Poil)pi + Pi{l)po' 



(43) 
(44) 



and the maximum of Qi is 

Qlmax ^1 



H{Po{l)po + Pi{l)pi,Po{l)pi + Pi{l)po) 
-H{po,pi) 



(45) 

Thus, it is found that Qimax is always positive what- 
ever the insertion position is. Combining Eq. (j45p and 
Eq. ([37)1 . one gets the best efficiency as 



1/2 
l/L 




(1) Bose Case 



(2) Fermi Case 



Figure 5: (Color Online) Work extracted W as the function of 
insertion position l/L and the working temperature Ti. (1) Bose 
case with A'^ = 4. (2) Fermi case with A'^ = 4. 



errors. Its initial state is 



PdO 



0.7000 |0) (0| 
0.0630 |2) (2| 
0.0081 14) (4| 



0.2100 |1) (1 
0.0189 |3) (3| 



(48) 



The temperature of the heat sink is fixed that T2 = 1. 
Denote the higher temperature and the insertion position 
as Ti and I respectively. Then the work extracted by 
the heat engine during a thermodynamic cycle can be 
numerically calculated for both the Bose case and the 
Fermi case. The results are shown in Fig. [5] As one can 
see, the patterns in Fig. [5] for different distributions are 
quite different. In fact, in the next section, it is found 
that the difference between the Bose case and the Fermi 
case is most significant in low-temperature regime. 



Vmax = 1 - „ (46) 

Imax 

= 1-^- (47) 
J-i 

Then it is concluded that for a single-particle SHE, one 
can always choose a good expansion position to let the 
SHE reach the Carnot cycle's efficiency even if the de- 
mon contains some errors. Recall that in the step 4 of the 
thermodynamic cycle, the removing process, one needs to 
thermalization the working substance firstly. However, 
for a single-particle SHE, it is found that the expansion 
position can always be chosen that the thermalization 
process is not needed. Thus, a single-particle SHE can al- 
ways be a reversible heat engine with two heat baths, and 
its efficiency reaches the Carnot cycle's efficiency. On the 
other hand, for a multi-particle SHE, the thermalization 
process is always needed in the removing process, which 
makes the heat engine irreversible. Thus, a multi-particle 
she's efficiency can not always reach Carnot cycle's ef- 
ficiency. 

At the end of this section, a concrete example is pre- 
sented to show the difference between the Bose system 
and the Fermi system. In this example, a four-particle 
system is considered as the working substance. The de- 
mon we chosen is a five-level quantum system with some 



V. LOW-TEMPERATURE BEHAVIOR WITH 
DIFFERENT QUANTUM STATISTICS 

In this section, the behavior of a multi-particle SHE 
when the working temperature Ti is very low is studied. 
In Sec. llVi a necessary condition for the heat engine to 
extract positive work is found that Ti > T2. Thus, if Ti 
considered is very small, the temperature of the heat sink 
will be T2 = 0. The demon in this case is free of errors 
because it is natural to get a pure initial state if the 
demon is re-initialized by a zero-temperature heat sink. 
The motivation to concern the demon without error is 
to see the difference between the Bose system and Fermi 
system more clearly without the effect of the demon's 
errors. As mentioned in the step 5 of Sec. IIIIl it is less 
likely for one to re-initialize such a demon by using a 
reversible scheme; however, one can use an asymptotic 
scheme to re-initialize the demon. 

Now, we describe the asymptotic scheme in details 
firstly. Similar to the re-initialization scheme in Sec. IIIH 
the asymptotic scheme also contains three sub-steps. The 
three sub-steps are illustrated in Fig. [HI In the first sub- 
step, the energies of the levels \i) for « are adjusted 
from Ai to a very small value A(7^ 0) adiabatically. Dur- 
ing this sub-step, there is no heat transferred. In the sec- 
ond sub-step, the demon is in contact with the heat sink 
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A, 



1^ 

Ai 



k_0 
(d) 



A 



step3 



1 = 2 Pr' = o 



A A 



0^. , 

Po = 1 



(c) 



Figure 6: (Color Online) Asymptotic erasure scheme for T2 
0. 



at temperature T2 = until it is in equilibrium with the 
heat sink. During this sub-step, there is no work applied 
and the finial state of the demon is a pure state for its 
temperature is zero when A ^ 0. In the third sub-step, 
the energies of the demon's levels are re-adjusted adia- 
batically to its initial value. Thus, these three sub-steps 
complete the re-initialization process of the demon. If we 
use this scheme, Qi in Eq. (1361) will not change and Q2 in 
Eq. ([57)1 will be revised to be the heat transferred in the 
second sub-step. Due to the irreversibility of the second 
sub-step, this scheme is not reversible. Nevertheless, if 
one lets A tend to 0, Q2 will tend to 0, which makes the 
scheme reversible. Thus, Eq. ([55]) and Eq. ([57]) will still 
hold if one uses this asymptotic scheme. In the following, 
this asymptotic scheme is used in which the relationships, 
Q2 = and = Qi, hold. 

Since the initial state of the demon is a pure state, 
which corresponds to the conditions T2 — and pj — 
6jQ, the work extracted can be given by the relationship 
W = Qi and dMD as 



N 



W = Ti 



(49) 



This is just the same as that in Ref. [6|] and it can be 
calculated by a different way as in Ref. [l^l- Then we 
will use Eq. (|49l) to study the Bose system and the Fermi 
system separately. 



First of all, we study the behavior of W/Ti when Ti 
0. According to Eq. (|i^ . one gets 



W 



N 

E 

i=0 



(50) 



Here, the expansion position li have been optimized as 
in Sec. IIVI Thus Pi{li) is always no less than Pi{l). For 
a deterministic distribution, {Pi{l)\Pi{l) — Sij}, W/Ti in 
Eq. (Isni) becomes to be \nPj{lj). Due to the optimiza- 



tion of Ij, we have 1 ^ Pji^j) ^ PjQ) = 1: or, Pj{lj) = 1. 
Thus, it is concluded that W/Ti is always when {Pi{l)} 
is a deterministic distribution. In other words, the nec- 
essary condition for a non-zero W/Ti is that the infor- 
mational entropy H{{Pi{l)}) of the distribution {Pi{l)} 
is non-zero. 

Using the joint entropy theorem [l3|, one gets 



s^s^s = H{{p,{i)}) +Y,m)s^ 



(51) 



where 534113 is the entropy of the working substance 
after the insertion and Si represents the working sub- 
stance's entropy in the condition that the number of 
the particles in the right compartment is i. Due to 
the non-negative property of entropy, one finds 6*8.1118 1^ 
H {{Pi{l)}) . Therefore if the entropy of the working sub- 
stance is zero, one will find that H{{Pi{l)}) — and 
W/Ti = 0. For a system with its ground state non- 
degenerate, its entropy always tends to zero when tem- 
perature tends to zero, which is referred to as the third 
law of thermodynamics. Thus, it is found that, if the in- 
sertion position I makes the ground state of the working 
substance non-degenerate, one will always have 



lim — 

Ti^O Ti 



0. 



(52) 



However, if the insertion position / makes the ground 
state of the working substance degenerate, there will be 
a chance that W/Ti tends to non-zero. Such a position 
I is called a degenerate point. Thus, the relationship 
between the low-temperature behavior of W/Ti and the 
degenerate point has been established. The similar result 
is found in Ref. [Tsj . 

Then, we look for the degenerate points for both the 
ideal Bose system and the ideal Fermi system. For the 
ideal Bose system, all particles condense in the ground 
state of the system when temperature is zero. Thus 
the only chance for the ground state to be degenerate is 
that the energy of the left compartment's single-particle 
ground state is the same as that of the right compart- 
ment, which happens when the widths of the two com- 
partments are equal to each other. Thereafter it is found 
that there is only one degenerate point for the ideal Bose 
system, i.e., I = L/2. 

For a Fermi system, the particles satisfy the Pauli's 
exclusion principle. When temperature is zero, all the 
energy levels under the Fermi surface are occupied. Thus 



11 
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0.2 0.4 0.6 0.1 

l/L 



0.2 0.4 0.6 0.8 

l/L 



0.2 0.4 0.6 0.8 

l/L 



Figure 7: (Color Online) W/Ti as the function of the insertion position l/L and the working temperature Ti. (1) A Bose 
system with N = 2. (2) A Bose system with iV = 3. (3) A Bose system with TV = 4. (4) A Fermi system with N ^ 2. (5) A 
Fermi system with TV = 3. (6) A Fermi system with TV = 4. 
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conditions of the degenerate points can be written as 



{hiT:y _ {h{N + i-i)TTy 



2mP 



2m{L - If 



(53) 



(1) 



(2) 



Therefore the positions of the Fermi system's degenerate 
points are 



Figure 8: (Color Online) Degenerate ground states for a Fermi 
gas. When the Fermi surface contains more than one single- 
particle states, there is a chance that the ground state of the 
Fermi gas is degenerate. For example, the states in (1) and 
(2) are two ground states of the system. 



if the «-th single-particle energy level of the left compart- 
ment and the j-th single-particle energy level of the left 
compartment are both on the Fermi surface, the ground 
state of the total gas may be degenerate, see Fig. [Sj The 
index i and j should satisfy the relationship that — 1 
equals to the number of the particles TV. Thus there are 
TV situations corresponding to i = 1,2,- •• ,TV. There- 
after it is found that there are TV degenerate points for 
an ideal Fermi system with TV particles. As the two com- 
partments are modeled as two one-dimensional infinite 
square wells with widths I and L — I respectively, one can 
workout the spectrum of them as in Eq. (|A1I) . Thus the 



I 

L 



N + 1 



for i^l,2,--- ,N. 



(54) 



So far we have revealed the qualitative difference of the 
work extracted by the engine between the Bose type and 
the Fermi type working substance at low temperature. 
That difference in the quantum statistical behavior be- 
tween the Bose and the Fermi working substance will 
disappear as the temperature Ti increase. The numer- 
ical simulation in Fig. [7] quantitively illustrates the sta- 
tistical behavior at finite temperature. When Ti tends 
to zero, the numerical calculations confirm the above re- 
sults. Similar result is also found in Ref. Il3li. 



VI. CONCLUSION 

In summary, the multi-particle SHE has been stud- 
ied in a complete quantum framework, where the re- 
initialization process of the demon is discussed in detail. 
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The errors are introduced in the demon so that one can 
use a finite-temperature heat sink to cool down the de- 
mon while the traditional demon must be cooled down by 
a zero-temperature heat sink. For a single-particle SHE, 
there is a scheme that the whole cycle of the SHE can 
be revisible. We illustrated that the efficiency of such a 
scheme is indeed the efficiency of Carnot cycle's and other 
schemes have lower efficiency. And it is proved that SHE 
can not exceed the limit of the SLoT. 

On the other hand, the properties of the SHE with 
Bose or Fermi particles are studied. In the high- 
temperature regime, these two distributions make no dif- 
ference. Both of them reduce to the result of Maxwell 
distribution. In contrast, the Bose system and the Fermi 
system differ from each other in the low-temperature 
regime. The origin of the difference is the difference of 
the degenerate points of these two distributions. The 
SHE whose insertion position is a degenerate point ex- 
tracts more work in a cycle. If the temperature is abso- 
lute zero, the ratio of the work extracted to the working 
temperature will show some uncontinuous properties. 

This work was supported by NSFC through grants 
10974209 and 10935010 and by the National 973 program 
(Grant No. 2006CB921205). 



Appendix A: Multi-Particle Partition Function 

In this appendix, the explicit expressions of the multi- 
particle partition functions for the ideal Bose gas and the 
ideal Fermi gas are presented. Let Zn{l, (3) be the canon- 
ical partition function for n particles, where /3 represents 
the inverse temperature of the system and I represents 
some parameters of the system, e.g., in our paper I means 
the width of the infinite square well. To present the ex- 
plicit expression of Zn{l, the explicit expression of the 
single-particle partition function Zi{l,/3) is given firstly. 
Then the relationship between the multi-particle parti- 
tion function and the single-particle partition function is 
shown. 

The single-particle partition function does not depend 
on which distribution the particle satisfies. Once one 
knows the single-particle spectrum of the system, one can 
write down its canonical partition function. For example, 
the single-particle spectrum in a one-dimensional infinite 
square well with width I is 



condition, one gets 



Ei{l) 



{hiirY 
2MP 



1,2,3, 



(Al) 



where M is the mass of the particle and i is the index 
of the energy levels. Then the single-particle partition 
function of this system is 



Zi(/,/3) = ^e- 



(A2) 



i=l 



When it comes to mult i- particle partition function, we 
first take the two-particle system as an example. In this 



7bose(fermi) 



-I3E„ 



-2/9 £;„ 



(A3) 

where the indices (3 and I have been omitted for simplic- 
ity. Here, the summation, J2'^^iexp{—2l3Em), can be 
viewed as a single-particle partition function with spec- 
trum {2Em}- Thus, it is found that g^n be 
decomposed into some single-particle partition functions. 
In general, all the partition functions of the ideal gas, 
^bose(fermi)^ Can be decomposed into some single-particle 
partition functions by this ways. The next task is to 
workout the coefficients in the decomposition. 

In the following discussion, the system considered has 
discrete spectrum and the index of the energy levels is 
denoted by i = 1,2,3,---. Let be the number of the 
particles in the i-th level. To be convenient, we denote 
exp{—l3Ei) as Xi for i — 1, 2, • ■ • and the multi-particle 
partition function becomes to be 



{ni}i=l 



(A4) 



Here, the summation condition is that " i = 

n, = 0, 1, 2, • • • for Bose case or rij = 0, 1 for Fermi case. 

It is not difficult to show that 

OO -J 

H{t) ^ = n 

n—O I 

oc 

E{t) 4 ^ z^™'r = n(i + ^it)- (A6) 

n=0 i 

In fact, ift = exp(— one will find that H{t) and E{t) 
arc nothing but the grand canonical partition functions 
for ideal Bose case and ideal Fermi case respectively. In 
mathematics, H{t) and E{t) are the generating functions 
for Z^°^'^ and Z^''""' respectively and they have the fol- 
lowing simple relationship 



H{t) 



1 



Ei-t) 



(A7) 



Next, the single-particle partition function with spectrum 
{jEm}, I]m=i cxp(-j/3£;„i) = Y^i^i^ is denoted as Pj 
for J = 1, 2, • • • . Its generating function is 



P{t) 4 E^:'*'=EE(^^*)' 



OO 



i j=l 

Thus, it is found that 

P{t) = 



Xit 



t dH{t) 



(A8) 



(A9) 
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i.e., 

By indefinite integration, one finds that 

oo ^ 

J2 -Pjt^ = lni/(i) + C, (All) 

j=i 

where C is the constant of integration and it is not diffi- 
cult to find that C = 0. So, there is the following impor- 
tant relationship 

oo ^ 

Compare the coefficient of in both the two sides of 
Ea. (IA12p and one can find that 

^r^ = E^p({^4), (A13) 

where {iq} — (ii,J2,''' i*™) and the summation condi- 
tion is "Ylia ~ ^ which corresponds a summation for 
all the Young diagrams with n boxes. The z{{ia}) is 
defined as 

z{{i^})=h\r' -m'^-^-iM-, (A14) 

and the P{{iot]) is the product of some single-particle 
partition functions as 

P{W})=Pl^ -P',^ (A15) 

I 



Thus we have decomposed the multi-particle partition 
function into some single-particle partition functions for 
Bose case. For Fermi case, using the relationship in 
Ea.(|X7|. one finds that 

E{i)^cMY.-^—P,tn- (A16) 

Also, compare the coefficient of in both the two sides 
of Eq. (|A16P and it is found that 



where the summation condition is the same as in the Bose 
case. 



Appendix B: The Upper Limit of Qi 

In this appendix, the upper limit of Qi in Eq. (j36p is 
estimated. To this end, a two-particle SHE is used as 
an example firstly, in which the demon considered is of 
three-level. The function is chosen as fi[j) = i + j 
(mod 3). Thus the heat absorbed from the heat resource 
is the one shown in Eq. (|39p . To estimate the upper limit 
of Qi, Eq. ([5^ should be re- written as 



Qi =TiH{{P,{l)}) +Fo(0PolnPo,o + A(0P2lnPi,o + P2(0PilnP2,o 
+Poil)pi lnPo,i + ^(OPolnPi.i +^2(0p2 111^^2,1 

-t-Po(/)p2 lnPo,2 + Pi{l)pi lnPi.2 + P2(0PolnP2,2, (Bl) 



where Pi j satisfy the constraints Pij — Pi{lj) for all 
i and j. Thus the original optimization problem be- 
comes to be an optimization problem with constraints. 
This optimization problem with constraints has a maxi- 
mum solution Qimax since Qi in Eq. (jBl[) has an upper 
bound TiH{{Pi{l)}) and the domain of {Pij} is closed. 
Substitute these constraints by the weaker constraints 
'YliiPij = 1 for all Jj then we get another optimiza- 
tion problem with constraints. After solving this prob- 
lem we will get another maximum Qimax- Because the 
constraints in the latter problem are weaker than the 
constraints in the former one, we have Qimax ^ Q'lmax- 
Thus Qiiuax is a upper limit of Qi. The calculation of 
Q'lmax can be performed in terms of Lagrange multipliers. 



Namely, one only needs to solve an optimization problem 
without constraint in which the target function is 

Ql{{P^.o}) + Ao E + ^1 E + ^2 E ^^2, (B2) 

i i i 

where Aq, Ai and A2 are three Lagrange multipliers. Solv- 
ing this constraint-free optimization problem, one finds 
that the maximum solution Q'lf^ax is reached when 
satisfies the following relationships 

Po,o : Pi,o ■■ P2,o = Po{l)po ■■ Pi{l)P2 ■■ P2{l)pi 

Po.l ■■ Pl,l ■■ P2,l = Poil)pi ■■ Plil)po ■■ P2{1)P2 (B3) 

Po,2 : Pi,2 : P2,2 = Po{l)P2 ■■ Pi{l)pi ■■ P2{l)pa- 
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For a general situation, this trick can also be used to 
find an upper limit of the heat absorbed from the heat 
resource Qi. To this end, one only needs to maximize Qi 

I 



= Ti 
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